Harnessing the GaAs quantum dot nuclear spin bath for quantum control 
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We theoretically demonstrate that nuclear spins can be harnessed to coherently control two- 
electron spin states in a double quantum dot. Hyperfine interactions lead to an avoided crossing 
between the spin singlet state and the m s = +1 triplet state, T+. We show that a coherent su- 
perposition of singlet and triplet states can be achieved using finite-time Landau-Zener-Stiickelberg 
interferometry. In this system the coherent rotation rate is set by the Zeeman energy, resulting in 
~ f nanosecond single spin rotations. We analyze the coherence of this spin qubit by considering the 
coupling to the nuclear spin bath and show that Tj ~ 16 ns, in good agreement with experimental 
data. Our analysis further demonstrates that efficient single qubit and two qubit control can be 
achieved using Landau-Zener-Stiickelberg interferometry. 



I. INTRODUCTION 

Considerable effort has been made in the past few 
years to implement qubits in nanoscale solid state struc- 
tures. One of the most promising candidates are spin 
qubits confined in electrostatically defined quantum dots 
(QDs) embedded in GaAs structures^ 2 -. A universal set 
of quantum gates has been demonstrated in GaAs dou- 
ble quantum dots (DQD) through the achievements of 
single-spin rotations and the two-spin exchange interac- 
tion that generates the VSWAP gate^— . Despite these 
advances, coherence times are limited by the hyperfine 
interaction, which couples the trapped electron spin in 
the quantum dot to the spin-| nuclei of the GaAs sub- 
strate. The resulting nuclear fields cause rapid electron 
spin dephasing, leading to an inhomogenous dephasing 
time T 2 * ~ 10 ns. As each electron spin is coupled to ap- 
proximately one million nuclei, the resulting behavior of 
the coupled electron-nuclear spin system is complicated 
and leads to rich dynamics that are sensitive to experi- 
mental parameters^. 

The hyperfine interaction has traditionally been 
viewed as a nuisance. However, a recent experiment 
demonstrates that generation of a controlled nuclear field 
gradient can be used to drive fast spin rotations^. The de- 
velopment of quantum control methods in semiconductor 
quantum dots that are based on nuclear spin interactions 
could lead to new paradigms for single spin control. 

We theoretically show that hyperfine interactions can 
be harnessed for quantum control in GaAs semiconductor 
quantum dots. In the presence of an external magnetic 
field B, which splits the triplet states, the hyperfine in- 
teraction results in an avoided crossing between the spin 
singlet S and spin triplet T + , which form the basis of a 
new type of spin qubit. Coherent quantum control of this 
qubit has already been experimentally achieved through 
Landau-Zener-Stiickelberg (LZS) interferometry^, a tech- 
nique previously used to coherently manipulate supercon- 
ducting qubits^^ and which is based on the interference 
due to repeated LZS tunneling events^—. The original 
LZS problem studies a two level system which exhibits an 
avoided crossing when an external control parameter is 
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Figure 1. (color online) (a) Energy diagram for the relevant 
states in the DQD as a function of e. The spin states for 
the implementation of the qubit are the hybridized singlet S 
and the triplet T+ . (b) A coherent superposition of the qubit 
states is generated by LZS interferometry. (c) The initialized 
S(2, 0) is swept through the avoided crossing by means of an 
applied linear gate voltage pulse ej —5- Ef. The final state 
is a coherent superposition of S and T+ generated by LZS 
tunneling. At et, the system evolves in the external magnetic 
field B for a time t w before a reverse gate voltage pulse brings 
the system back to e\ where a QPC measurement is performed 
to determine the singlet state return probability, Ps. 



changed. If the control parameter is time dependent, the 
system can be brought from an initial state through the 
avoided crossing. This passage may result in a change of 
populations and relative phase of the states. In the LZS 
version, the system is driven from t\ = — oo — > t{ = oo 
and the difference in energy between the two states is a 
linear function of time, AE(t) = at, which leads to the 
well-known result for the non-adiabatic transition prob- 
ability P LZS = exp(-2^|(S| J ff int |T + )| 2 /aft). 

For the DQD system, the usual infinite-time asymp- 
totic theory describing LZS interferometry cannot be 
used. The avoided crossing originates from the hyper- 
fine interaction between the electronic spins and the nu- 
clear spins whose fluctuations result in a poorly defined 
crossing position. The phenomena observed in the exper- 
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iments can only be properly described using finite-time 
LZS theoryi 5 .. 

To prove that our formalism describes correctly the 
coherent manipulation of the S-T + based qubit, we will 
compare it to the experimentally measurable quantity 
Ps, the singlet return probability. We then show how 
single qubit operations can be engineered either by the 
Euler angle method for rotations or by only using LZS 
interferometry. Finally, we demonstrate that a two qubit 
gate can be achieved by capacitively coupling two S-T + 
qubits. In contrast to the S-T qubit, where the rotation 
rate is set by a charge-noise-susceptible exchange energy, 
the rotation rate in the S-T + qubit is set by the Zeeman 
energy and approaches 1 ns for modest magnetic fields of 
100 mT. 



II. MODEL 

The spin preserving Hamiltonian, 
Ho = ^e js n js +u^n n n ji +T^(^c[ s c 2s +h.cj, (1) 

js j s 

with Ej S — Ej + g* hbBs/2 and rij s = c^ s Cj s describes 
the coupling between two electrons in a DQD in a mag- 
netic field B. g* denotes the effective Lande g-factor 
(—0.44 for GaAs), /zb the Bohr magneton, the j = 1,2 
and s =t,4_= ±1 label the dot number and spin. The 
first term is the single-particle energy of the confined elec- 
trons, the second accounts for the Coulomb energy u of 
two electrons on the same QD, and the last for tunneling 
with strength r between the dots. 

The diagonalization of the first two terms of Eq. ([1]) 
leads to the relevant charge states of a DQD: the sin- 
glets S(0, 2), S(2, 0), S(l, 1) and triplets T±, (l, 1) where 
(I, r) denotes the charge configuration of the dots (see 
Fig. [Ha)). The other states can be neglected as they have 
energies much higher than those considered here. The de- 
generacy of the singlets S(2,0) and S(l, 1) at e = ±u is 
lifted by the inter-dot tunneling, resulting in a splitting of 
y/2r. The hyperfine interaction Hup = Si • hi + S2 • h\2 
between the electron spins Si and the nuclear spins \\ 
opens a splitting Ahf at the degeneracy point e c of the 
singlet S and the triplet T+. Here, hj = ~£l=\ A i ! i is 
the Overhauser (effective nuclear) field operator. The 
sum runs over the n(i) nuclear spins in dot i and A\ — 

VikVo \^i{ T k)\ is the hyperfine coupling constant with 
the fc-th nucleus in dot i, with Vt'i(rfc) the electron wave 
function, vq the volume of the unit cell and Vik the hy- 
perfine coupling strength. From now on, since we assume 
symmetric dots, we have n(l) = n(2) — n. Introducing 
Sf = Sf ± iSf and hf = hf ± ih%, we write 

#hf = \ J2 (2 W + S+K + Sr h+) . (2) 

i 

To determine which spin states are relevant for our 
theory, we consider Plzs from the asymptotic LZS model 



and the result of Ref. [16| for multiple level crossings to es- 
timate the order of magnitude of the transition probabil- 
ities. The initialization of the system is done by prepar- 
ing a singlet S(2, 0) (ej > e c ), then e is swept to achieve 
Ef < e c . During this operation, the system goes through 
three avoided crossings (cf Fig. [T]). To estimate the ma- 
trix element | (S|#hf|T+) | 2 = |(S|P HF |T_)| 2 = |A HF | 2 , 
we use the experimentally found value of Ahf = 60 neV 
from Ref. 0. The matrix element entering Plzs at the 
avoided crossing between S(2, 0) and the excited singlet 
state S'(l, 1) is given by |(S(2, 0)|£T |S'(1, 1))| 2 = 2r 2 , 
with t = 5 fieV. The order of magnitude of a is taken 
between 10~ 3 - 10~ 2 meV/ns. We find P T+ ~ 0.97 - 0.99 
and Ps'(i, i); -Pt_ *C 10~ 8 . These results show that pop- 
ulation of the excited singlet and T_ level are negligible. 
To can also be neglected because it does not cross with 
any other level, it splits from S(l, 1) due to the exchange 
coupling^ 2 -. 

Near the S-T + crossing, the dynamics can be restricted 
to the Hilbert space spanned by T + (l,l), S(l,l), and 
S(2, 0) and described by 

(g*HnB \ 
#S,t + ^ V2r (3) 
\ V2t u-e) 

where we can neglect an additive term oc el, with 
c = E\ — £2 the detuning of the dots. According 
to our previous estimate, we can reduce Eq. Q to 
a 2 x 2 Hamiltonian which only takes into account 
the lowest hybridized singlet |S) = c(e)|S(l, 1)) + 
s/l - c(e-) 2 |S(2, 0)) and the triplet |T+) = |T+(1, 1)), 
where c(e) = (— u + e — 77) /-\/8t 2 + (u — e + tj) 2 , with 
i] = ^/8t 2 + (u — e) 2 . The energy associated with the 
lowest hybridized singlet is Es(e) — (u — e — rf) /2 and 
the energy of the triplet is = g* ^l-qB. The Hamilto- 
nian describing the dynamics of the lowest energy states 
in the vicinity of S-T + can therefore be written as 

H (e) = E S (e) I S) (S I + E T+ \ T+ ) (T+ 1 . (4) 

Another relevant quantity derived from Eq. ([3]) is the 
degeneracy position e c of S and T + , given by the funnel- 
shaped function 

e c (u, B)= u + 2r 2 /gfi B B - gfi B B. (5) 

III. SEMI-CLASSICAL THEORY 

We model the Overhauser field classically, such that it 
acts on the electron spin as a magnetic field B ni , with 
hi = g*/iBB n! i and its physical properties are given by a 
statistical distribution that reflects the quantum fluctua- 
tions of the nuclear ensemble. At typical operating tem- 
peratures and external magnetic fields ksT 3> gtq/J,^B, 
where and fi^ are the nuclear g factor and magne- 
ton. In this limit we can assume the nuclei to be com- 
pletely unpolarized, resulting in a Gaussian distribution 
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of nuclear fields^ p(S n>i ) = (l/v / 2^a)e- B "^ /2 ' T , with 
cr = A/ g* HBy/n and A ss 90 //eV. The effective Hamil- 
tonian describing the qubit dynamics around the S-T + 
avoided crossing is given by 



H cS (e) = H (e) + \g* m £ (S^B^ + 5r5+, 



(6) 



We include the nuclear Zeeman spliting in the T+ en- 
ergy, P T+ = g*p B (B + B*), where B* = B* x + B* a 
and B^ i = B* i ± iB^ f . The classical approximation 
of the Overhauser fields is possible because the nuclear 
state changes only slightly after a single sweep^. To ob- 
tain an analytical expression of the LZS propagator (see 
Appendix IX)) . we linearize the difference in energy AE(t) 
around t — 0— by assuming e(t) — jt + £ c , we find 



2r 2 + ( 9f x B (B + B%)} 



r7- 



(7) 



Here 7 is the rate at which the external voltage gates 
are ramped. Control over a can therefore be achieved by 
modifying 7. 

To test our model we first compute the singlet return 
probability Ps , an experimentally observable quantity, as 
a function of the final detuning S{ and waiting time t w , 



Ps = f I] dB n , fcP (B n , fc )|(S|[/(B n , 1 , B n , 2 )|S)| 2 , (8) 

J fe=l,2 

with 

C/(B n ,i,B n , 2 ) = P b (B n ,i,B n . 2 )P w (B*)[/ f (B n ,i,B n , 2 ), 

(9) 

where Ub,{ = T exp[— i J. f dtH c g(e(t))/h] are the 
backward and forward LZS propagators and U w ~ 
Texp[— i J W dt/fo(ef)/ft] describes the evolution of the 
system during the waiting time £ w at the final detuning 
position £f with \ef — £ c | > Ahf- We evaluate Ps by nu- 
merical sampling. Instead of estimating a from n and A 
of the QDs, here we use the experimentally determined^ 
A HF = 60neV to derive a ~ 1.67 mT. We have (A^ F ) = 
((SB-) 2 

(t+ip; 

and t w for 7 = 0.015 meV/ns, B = 100 mT, u = 4meV, 
and r = 5 /ieV. We use a square pulse with a ramping 
time |tf — t}| fixed to 1.5 ns and the initial detuning e; is 
varied to reach different values of £f. 

We identify coherent oscillations as a function of f w 
and £f. From a best fit, we obtain the decoherence 
time T 2 * = (16.0 ± 0.4) ns, which agrees well with ex- 
periment. The decoherence is mainly due to the fluctu- 
ations of B* . The period of the temporal oscillations is 
T = h/ \Es{e t ) - E T+ \ - 4.3ns for e f = 3.97 meV (see 
Fig. E[b)). For a fixed B, a shorter period can be ob- 
tained for smaller ef, the fastest oscillations being defined 
by the Zeeman energy. To further decrease the period the 



- 4a 2 , with SB- = g*m(B- 2 - B~ x )/2y/2 = 
hf|S). In Fig. [5] we show Ps as a function of £f 
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Figure 2. (color online) Theoretical results, (a) The sin- 
glet return probability Ps from the semi- classical model as 
a function of the waiting time t„ and final detuning Ef. We 
find nanosecond oscillation periods and the dephasing time 
T 2 * ~ 16 ns, in good agreement with experiment. We used 
B = 100 mT, 7 = 0.015 meV/ns, u = 4meV and r = 5/ieV. 
(b) Ps (blue) as a function of t w for et = 3.97 meV (hori- 
zontal line in (a)). We plot (orange) Csin(uit) exp( — (t/T^) 2 ) 
where w = |P S (ef) - E T+ \ /h ~ 2vr ■ 0.23 GHz, C = 0.95 and 
T2 = (16.0 ± 0.4) ns is extracted from a best fit. 



external B field can be increased and hence the qubit 
manipulation could be done in a time scale of 100 ps for 
B ~ 1.6 T, which would allow ^160 coherent operations 
within T 2 . In the exchange gate demonstrated in Ref. [3, 
dJ/de increases with e, which results in faster dephasing 
for faster rotations. In contrast, here the rotation rate is 
set by the Zeeman energy, which is independent of e far 
from the S-T + avoided crossing. As a result, the coherent 
oscillation frequency can be increased without making 
the qubit more susceptible to gate voltage fluctuations 
by simply increasing B. Far from the avoided crossing 
the level spacing is independent of detuning, similar to 
the "sweet spot" in superconducting qubits 2 ^. 

The model predicts coherent oscillations in Ps for 
Ef > e c , i.e. in the case where the qubit has not passed 
the avoided crossing. It can be explained within finite- 
time LZS theory, but not with the conventional asymp- 
totic Landau-Zener formula. In other words, even if 

if < 0, we have \Ujf s (tf, ti)\ 2 > 0, which illustrates the 
non-adiabatic character of the problem. For the pulse 
conditions used in Ref. 8, oscillations are not observed 
for £f > e c , most likely due to charge dephasing. The co- 
herence time of an admixture of (1, 1) and (2, 0) charge 
state has been measured to be ~ 1 ns for GaAs QDs2i and 
sets the time scale at which the system must be driven 
to observe oscillations for ei > e c . A finite-time effect in 
agreement with the experimental data is the dependence 
of the oscillation amplitude Ps on the pulse length. Fi- 
nally, we show in Fig.[3]a comparison between experiment 
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Figure 3. (color online) Comparison between experimental 
(purple, open circles) and theoretical values (blue, filled cir- 
cles) of Ps for B = 45 mT. Theoretical values are obtained 
by finding the detuning et = 3.78 meV for which the theoret- 
ical and experimental oscillation periods match. A suitable 
7 = 0.12 meV/ns is then chosen such that gj > e c . The exper- 
iments used passive filtering of square pulses to reduce Plzs 
and thereby increase the visibility of the oscillations. As a 
result, oscillations in the experimental data are delayed to 
longer times due to finite rise time effects. The theory points 
are obtained with a perfect square pulse. 

and theory for B = 45 mT. The experimental data were 
obtained from the setup used in Ref. [H. The experiments 
used passive filtering of square pulses to reduce Plzs an d 
thereby increase the visibility of the oscillations. As a 
result, oscillations in the experimental data are delayed 
to longer times due to finite rise time effects. 

IV. ARBITRARY SINGLE QUBIT ROTATIONS 

The passage through the avoided crossing can be inter- 
preted as a rotation (see Appendix [Cjl . U^zsiVi t{, t{) — 
c -ier-n0/2 an an gi e q around the axis n, see Fig. |U 
Here 9 = 9(rj, tf, ti) and n = n(rj, tf, t{) is a unit vector 
where rj — |(T + |P/hf|S)| j \fah is the dimensionless cou- 
pling strength. Since n and 9 are functions of the same 
experimental parameters £;(—>• £;), tf(— > £f), anda(— > 7), 
it is not straight forward to find them simultaneously in 
order to build a given single qubit rotation. For instance, 
fixing two parameters and tuning the third one will si- 
multaneously change h and 9 limiting the achievable ro- 
tations angles. Nevertheless, the situation is not hope- 
less and several composite methods can be engineered to 
achieve any rotation. We present here three methods, 
each of them having their own advantages. 

Since rotations by an angle <p around the z-axis are 
available by letting the qubit evolve in an external B 
field, we would like n to be in the xy-plane in order to 
build any rotation by the Euler angle method. Below we 
show that this is possible if for example (see Appendix 
[Cjl the propagation times are equal, — tj = tf = £lzs- 
However, a 7r-rotation from |S) to |T + ) would take an ex- 
ponentially long time with a single LZS transition, since 
it corresponds to a fully adiabatic transition. However, 



Figure 4. (color online) Bloch sphere representation of LZS 
transitions as rotations, (a) When the propagation times are 
equal, —ti = tf, the rotation axis lays in the a;j/-plane. In 
addition with the rotations generated around the z-axis by 
letting the qubit evolve in the external magnetic field, any 
rotation can be achieved by the Euler angle method, (b) If 
the propagation times are different the rotations axis may not 
lay (see Appedix [C|) in the xy plane. In this case, a pure LZS 
interferometry technique can be used to generate any rotation. 



this problem can be circumvented by sequentially apply- 
ing several LZS transitions. A 7r-rotation from |S) to |T+) 
can be achieved in ~ 0.1 — 1.0 ns for two consecutive and 
identical LZS transitions. 

The single qubit gates can also be implemented with 
a pure LZS interferometry technique, similar to the one 
used to control superconducting qubits. This method re- 
quires sequential driving of the qubit through the avoided 
crossing. The different passages result in a series of LZS 
transitions each of them corresponding to a rotation of 
the qubit. By tuning 7 and choosing — tj 7^ if with differ- 
ent ratios for |fi|/|tf |, any qubit rotations can be achieved 
within a nanosecond. 

Since finite-time effects are present in the system, we 
can think about a control method where the qubit is op- 
erated on the (1,1) charge configuration side ("sweet re- 
gion"). This requires the preparation of a T + (l, 1) state 7 
and pulses with rise times shorter than 1 ns which do 
not drive the system through the avoided crossing unless 
a measurement is required. The qubit manipulation is 
achieved through finite-time LZS interferometry, where 
tuning 7 and the propagation time t- h f allows to achieve 
any desirable angle. For all those methods, an arbitrary 
single qubit rotation can be expressed as a series of a 
forward sweep -wait -backward sweep operators, 

v{B, V ) = l[u b (o®)u,(<p®)u f (e®), (10) 

i=i 

which reduces to Eq. ([5]) for I = 1. The proposed methods 
require a maximum of I = 3. It is important to notice 
that the rotation axis and the final measurable angles will 
not be n, 9 = Y^i^ii an d V — J^i^^ but rather (n) x , 
(0) x , and (<p) x , where the brackets denote the averaging 
over the nuclear spin bath. A similar scheme with / = 1 
has been proposed for the S - Tq qubit^. 
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Figure 5. (color online) A conditional gate can be imple- 
mented by capacitively coupling electrons trapped in quan- 
tum dots belonging to different qubits. The crossing position 
£ c has different values whenever the charge state of the con- 
trol qubit is (0 2) (a) or (1, 1) (b). The later case results in 
£f <C £c which suppresses any LZS transition. 



V. TWO QUBIT GATE 

To complete the set of quantum gates, a two-qubit 
operation such as CNOT is required. We consider the 
Hamiltonian 



H. 



(i) 



H. 



(2) 



(11) 



(i) 

where H^J is the single qubit Hamiltonian (JoJ) and H u 
U ^2j—2 3 n jf n jl describes the capacitive coupling^ 
between two adjacent QDs belonging to different S — T + 
qubits. Tunneling between the dots of the different qubits 
can be suppressed by an appropriate gate voltage. If the 
control qubit, qubit-1, is in a S(2, 0) state, H int = and 
the dynamics of the target qubit, qubit-2, is reduced to 
the case of a single qubit, see Fig. [5] When the control 
qubit is in a (1, 1) charge configuration, the target qubit 
is influenced by the interdot Coulomb interaction. In this 
case, the dynamics of the target qubit can be described 
by Eq. ([5]) by replacing u — > u + u. In particular, this 
affects the position of the avoided crossing. For a system 
with two DQD separated by a distance d ~ 2R — 10R, 
where R is the approximate radius of one QD, the in- 
tradot Coulomb interaction u ~ e 2 /R is comparable to 
the interdot Coulomb interaction u ~ e 2 /d resulting in 
\e c (u) — e c (u + u)\ > 1 meV. From the previous discus- 
sion, we know that a 7r-rotation is possible within ~ 1 ns 
if e c = e c (u). In the case where the avoided crossing is 
at e c (u + u) the same LZS sequence will leave the tar- 
get qubit unchanged, even within a finite-time theory 
since the separation between the two avoided crossings is 
> 1 meV. Therefore, we estimate the CNOT gate time 
to be ~ 1 — 3 ns. 

Let us consider the case where the control qubit is in 
the S(2, 0) state, which is the logical |1) of the qubit, such 
that the target and control qubit are not capacitively 
coupled, u = 4meV. For 7 = 0.3 B = 100 mT, 
£i = 3.98597 meV, and e f = 3.97991 meV we find for 
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I = 6 and tp® = in Eq. HH 



-0.314 0.945 - 0.096i\ 

0.945- 0.096i -0.314 ' " ^ ' 



This example shows the almost perfect realization of a 
conditional ia y operation which corresponds to a CNOT 
gate up to single-qubit gates. 

To show that this method produces a CNOT gate, 
we consider the case where the control qubit is in the 
T + (l, 1) state, which is the logical |0) of the qubit. We 
estimate a lower bound for the strength of the capacitive 
coupling between the qubits to be u + u = 5 meV (see 
above). In this case, the target qubit evolution takes the 
form 



P(°>(0, <p) 



0.971 + 0.002i -0.206 
0.206 - 0.124i 0.971- 



-0.124i\ , , 
0.002i 1 ' { ' 



which is close to 1 and demonstrates the possibility of 
generating a CNOT gate with the proposed method. 

Notice that our choice for s- h f corresponds to a prop- 
agation time tLzs — 0.01 ns such that the total gate 
time to achieve controlled-i<7 y is 0.24 ns. The fidelity is 

T = |TrX>(°) +Tr [V^{ia y Y]\ 2 /16 ~ 0.918. A more 
accurate CNOT gate can be engineered by fine tuning 
the parameters entering the LZS propagator. 



VI. CONCLUSIONS 

We have demonstrated that coherent control of the S- 
T_|_ qubit can be achieved using LZS interferometry. Hy- 
pcrfine interactions lead to an avoided crossing between 
S and T + states, which allows for efficient quantum con- 
trol. Moreover, we predict that in the limit of fast rise- 
time pulses coherent oscillations in Ps should be observed 
even without going through the avoided crossing. This 
phenomenon is a finite time effect which we have theoret- 
ically described using the general finite-time LZS theory 
and it can be used to operate the qubit in the (1, 1) 
charge configuration side ("sweet region"). 

Our scheme can be extended to DQD in materials with 
few nuclear spins (graphene, CNT, Si). In such cases, 
the avoided crossing between the qubit states can be 
achieved by engineering a DQD in the presence of micro- 
magnets which provide the in-plane gradient magnetic 
field for the realization of the LZS based gates 25 . The 
qubit will moreover benefit from the lack of the inho- 
mogeneous broadening due to the Overhauser fields and 
exhibit an extended T 2 *- l n GaAs DQDs the method 
proposed in Ref. could be used to extend T 2 * without 
cancelling the gradient field. Other schemes to polarize 
the nuclear spins^S or reduce their fluctuation o 27 ' 28 also 
exist. 
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Appendix A: The Landau-Zener-Stiickelberg 
finite-time propagator 

In this appendix we follow the work of Vitanov and 
Garraway^ and consider, without loss of generality, a 
two-level system whose eigenenergies E\ and E 2 are time 
dependent, E\ — E\(t) 1 E 2 = E 2 (t) and their difference 
is a linear function of time 2A(t) = E 2 (i) — E\(t) = at. 
Furthermore, we assume the levels to be coupled with 
strength A. The matrix representation of the system's 
Hamiltonian is given by 



#LZSW - ( A A(t) 



(Al) 



The time evolution of such a system is described by 
the time-dependent Schrodinger equation 



ih-m)) = H hZS (t)\m) 



(A2) 



with \tp(t)) = ci(t)|l) + c 2 (t)|2). After substitution of 
(|A1[) into Eq. (| A2|) . a coupled system of first order ordi- 
nary differential equations is obtained 



ifici(t) = -A(t)ci(i) + Xc 2 (t), 
iHc 2 (t) = Aci(t) + A(t)c 2 (t). 



(A3) 
(A4) 



By deriving Eq. (| A3|) with respect to time and sub- 
stituting Eqs. (|A3|) and ()A4[) into the newly obtained 
ordinary second order differential equation, we obtain 



c x (t) 



7?) 



(A5) 



It is convenient to introduce dimensionless parameters 
before solving Eq. (|A5|) . here we introduce the dimen- 
sionless time r = y/jrt which we substitute in Eq. 
to obtain 



dr 2 



Cl (T) + (-i + 77 2 +T 2 ) Cl (r) -0, (A6) 



-^=p is the dimensionless coupling strength. 



where r\ 

The solution of Eq. (TAH is 

ci(i) = KiD^ ^e^rj + K 2 D in 2 (y2e~* 



4 T 



(A7) 



where D v {z) are parabolic cylinder functions, which solve 
the Weber equation^ 



(A8) 



and can be obtained from Eq. (|A6j) by writing the ex- 
pression in brackets as — 2i(i77 2 /2 + l/2 + ir 2 /2) and sub- 
stituting t — > 2~ 1 / 2 exp(i7r/4)z. 

c 2 (t) is obtained by inserting Eq. (|A7I) into Eq. (|A3| 
and using the property 



d / £_ \ 

— [e 4 D v (z)J = ve * D v -i(z). 



(A9) 



One finds 



c 2 (t) = -^=e " -KxD^_ x (\/2e ^ t 



+K 2 D ij;l i (y2i 



e 4 r 



(A10) 



To find the constants n\ and k 2 , we consider initial 
conditions given by c\(t\) and c 2 (ti) and the Wronskian 
relation 

W {£>„(«), D v (-z)} := D v (z)^D v {-z) - D v (-z)^D„(z) 



'2tt 



r(-iz) 



(All) 



We solve the system of equation given by 
Eqs. (|A7|) and (IA10|) for ki and k 2 using the Wronskian 
property (|A11|) . we find 



- [^_ 1 (V2e^r i )c 1 (r i ) 



V2 



«2 



e~ D 



V2 ,„ 



^ ^\/2e 3 ^Ti^ c 2 (Ti) 

» i2 2_ 1 (v^e"^^ ci(Ti) 
^2 ^V2e~T Ti ^ C2 ( T; ) 



(A12) 



(A13) 



Substituting Eqs. (|A12[) and (|A13I) into 

Eqs. (|A7|) and (|A10|) and having in mind that we 

are looking for the evolution operator E/(tf, i;) giving 
the final state knowing the initial one 



|V(tf)) = Ufa, *0|V(*i)> 

we finally find the LZS propagator 



^LZ S (if, ti) 



Uu(t{, tj) U 12 (t{, tj) 
M 2 l(tf, ti) U 22 (t{, ti) 



(A14) 



(A15) 
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with 

un(t{, ti) = u* 22 {t u ti) 



where we have defined 



r(i-if 



/27T 



(V^e^Tf) %^ (VSe - ^ 



+ 77 x/ 1 " 2 + ?? 2 



and 



(B5) 



(A16) 



and 



cos 



"12 (*f, *i) = ~«2l(*f) *i) = 

(y^e-Trf'j (V^e^Ti 



\ 



a/t 2 + r; 2 



r(i-f 



e 4 



sm tf = 



\ 



(B6) 



2 ( X v^T^ 



(A17) 



Using the above expressions and writing n = e 17r |ri| to 
fulfill the condition for the expansion we find, 



In the original LZS problem t — is defined at the energy 
levels crossing. A situation where ti < and tf > corre- 
sponds to drive the system through the avoided crossing. 
The case tj < and tf < corresponds to stop the sys- 
tem before it goes through the avoided crossing. Finally, 
ti > and if > corresponds to a system which is ini- 
tially prepared after the avoided crossing. 



Appendix B: Asymptotic expansion of the parabolic 
cylinder functions 

The expression of the LZS propagator can be expressed 
with simpler functions when the argument t> 1 and the 
parameter r\ 3> 1, in this case the parabolic cylinder func- 
tions can be expanded asymptotically—. The necessary 
asymptotic forms to expand Eqs. (|A16I) and (|A17[) are 



«ll(*f, U) = «22(*f) *i) 



x/T 



e 7n ? sm Of cos t^e 



-i(€f-Ki+ar g r(l-ii- 



iUf+6+argr ll-SL. 



e 2 



r£ ( 

2 ( cos Of 



sm tf ; cos tff e 
cos0ie i(& ~ Si) - shift sin6>ie _i(&_Si) 



and 



(B7) 



Uu(tf, ti) = -«2l(*f) f i) - 

y/l e-^ (^sin^sin0 i e- i ( ?f+ « i+arsr ( 1 -*) + ^) 
-co S e f cosfte 1 M 1+argr ( 1 -^) + ^) 



L> ±i „2 fy^c^r) ~cos6>e^ ±i€ , (Bl) + e"^ (cos (9; sin 6» f e" i(?f 



cos tff sin fte 



C± i2i _ 1 (v^e^r) ~ — sinftr£ ±i <*+*> , (B2) 



,„ 2 fV2« 

2 \ 



ni-f 



sm 



(B3) 



(B8) 

where cos ft f, sin ft . f, and & f are respectively given by 
Eq. (JB6J) and Eq. flU) for r = e i7r | ti | , r f . 

We noticed that the asymptotic expansions (|B1[) . (IB2|) . 
(|B3|) . and (|B4[) are valid for the weaker condition r + r; S> 
1, as already reported in Ref. UH for the expansion of 
|un| 2 and \u 12 \ 2 . 



Appendix C: The LZS propagator as a rotation 

In quantum mechanics the rotation operator T>{n, •&) 
by an angle •& around an axis h of a two-level system has 
the representation 



D 



sin to 8 +H5 



V(n, tf) = e iflc 



foZ 2 (B4) 

■ cos 0e 8 4 , 



(l-f) 



cos 5 — m z sm ^ ( — m x — n y ) sm ~ 
(— iria: + n y ) sin $ cos ^ + in z sin 



(CI) 
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Figure 6. (color online) (a) Cosine of the rotation angle and 
components of the rotation axis (b), (c), and (d) as a function 
of Tf for the dimensionless parameters n = 10 and r\ = 3. 



Identifying Eqs. (jCip and (IA15[) with u^a given by 
Eqs. (|B7|) and (|B8|) we can express the rotation angle 
•d and the rotation axis h as functions of the LZS propa- 
gator parameters t\, T{, rj. We have 



cos 



r- 
77 = V 1 



7r?)2 cos £ [sin 6{ cos 6*; + cos 9f sin ( 
c 2 cos <r [cos 0f cos 6*i — sin 0f sin 9\\ , 



where 



(C2) 



and 



q = & - a- 



7T 

4' 



(C3) 



(C4) 



The components of the rotation axis are given by (see 
Fig.© 



- r -£ ( \/\ — e 7 "' 2 sin £ [sin 0f sin 9\ + cos #f cos 9[] 



sin - 



-e 2 sin g [cos #i sin 9i — cos #f sin 9i] I , 



(C5) 



— ~ a - f \/l — c 7r '? 2 cos C [— sin 0f sin 0; + cos 0f cos6*i] 
sin ¥ V 



+e 2 cos p [— cos 0; sin 0f — cos 9f sin 6?i] I , 



(C6) 

- — j ( y/l — c~ m ' 2 sin C, [sin 9{ cos 6*; — cos 9f sin 0i] 

2 

+e 2~ sin [— cos 0f cos 0; + sin 9f sin 0;]^ . 

(C7) 



smf 
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